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We argue that strong final state rescattering among the 
secondary particles created in relativistic heavy ion collisions 
is essential to understand the measured Bose-Einstein cor- 
relations. The recently suggested "random walk models" 
which contain only initial state scattering are unable to re- 
produce the measured magnitude and iG_-dependence of R± 
in Pb+Pb collisions and the increase of Ri with increasing 
size of the collision system. 
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An important aspect of understanding ultrarelativistic 
heavy ion collisions is the clear identification of genuine 
collective nuclear effects which cannot be explained in 
terms of a simple superposition of nucleon-nucleon colli- 
sions. An integral part of our task to look for new physics 
must therefore be the careful construction of models for 
nucleus-nucleus collisions ("A-B collisions") based on a 
superposition of individual nucleon-nucleon (N-N) colli- 
sions in order to establish where they fail. 

Recently there have been some renewed attempts to 
construct such models. In the "Random Walk Model" 
(RWM) of [Q the single-particle transverse mass spectra 
measured in A-B collision have been calculated by ex- 
trapolating those from pA collisions. The LEXUS model 
of j| goes even further, by simulating A-B collisions as 
a simple folding of independent N-N collisions. 

Here we demonstrate that, while these models have 
had some success in describing measured single-particle 
spectra, they fail to reproduce crucial features of the ob- 
served two-particle Bose-Einstein correlations. This is 
shown to be due to their lack of final state rescattering. 

The idea behind the RWM was to provide an alter- 
native interpretation of the measured single-particle m± 
spectra, opposite in spirit to the popular hydrodynamical 
parametrizations || . The latter are based on the assump- 
tion of a locally thermalized hadron resonance gas which 
undergoes longitudinal and transverse hydrodynamic ex- 
pansion. The transverse expansion is interpreted as a 
genuine nuclear collective effect with no analogue in N-N 
collisions; it is identified through the characteristic flat- 
tening it causes in the transverse mass spectra, especially 
at small M± < 2mp where the inverse slope parameter is 
found to increase linearly with the rest mass mo of the 
produced hadrons |||| . 

The RWM, on the other hand, starts from the obser- 
vation that such a flattening, relative to N-N collisions, 



happens even in pA collisions where hydrodynamic trans- 
verse expansion is not expected to occur. In the RWM 
an incident nucleon undergoes multiple collisions in the 
nuclear target, leading to a random walk pattern in the 
transverse momentum plane. When it collides inelasti- 
cally it creates a little "fireball" , identical to those formed 
in elementary inelastic pp collisions, but moving with a 
transverse rapidity p which is Gaussian distributed. In 
A-B collisions the same mechanism works for both pro- 
jectile and target nucleons. The width of the p distri- 
bution is fixed in pA collisions and then extrapolated to 
A-B collisions using geometric considerations Q. Reso- 
nance decays are neglected, limiting the applicability of 
the model to sufficiently large transverse momenta. 

In LEXUS the same idea is implemented more di- 
rectly in terms of a folding algorithm for N-N collisions 
which, in contrast to RWM, also includes longitudinal 
momentum degradation in multiple collisions. However, 
the latter is not coupled to the broadening of the re- 
distributions which is parametrized in terms of an effec- 
tive temperature (inverse slope) which again follows a 
simple random walk rule. 

In both models the transverse broadening parameter is 
taken as independent of rapidity and not correlated with 
the position of the collision point (fireball) within the 
reaction zone. Secondaries do not interact, and (except 
for energy degradation effects) all primary collisions are 
alike. For this reason neither model can explain the dif- 
ferent chemical composition in pp and A-B collisions (see 
iQl); they both make statements only about the shape of 
the momentum distributions. 

Existing tests of the models are still very superficial. In 
Q the LEXUS output was compared to pion and proton 
rapidity spectra and midrapidity transverse mass spectra 
from 200 A GeV S+S collisions at the SPS. In Refs. §§ 
recent data on single-particle m^-spectra from 160 A 
GeV Pb+Pb collisions at the SPS, taken by the NA44 
H and NA49 Q| collaborations, were compared to the 
RWM and to hydrodynamical parametrizations. Both 
types of models give a reasonable description of the data; 
discrepancies occur mostly at low p±, where resonance 
decays distort the pion spectra and hydrodynamical flow 
effects show up most strongly for the heavier hadrons. In 
particular, the RWM has difficulties || to reproduce the 
strong rest mass dependence of the transverse slope pa- 
rameters published by NA44 H which have been quoted 
as evidence for collective transverse flow It is clear 

from these papers, however, that without a systematic 



1 



study of the mj_-spectra as a function of the collision 
system and of rapidity, carefully including resonance de- 
cays, single particle spectra alone may not lead to definite 
conclusions about the validity of the RWM approach. In 
the present note we therefore examine the discriminating 
power of two-particle Bose-Einstein correlations. 

For ease of language our discussion will be based on 
the RWM, with comments on LEXUS to follow below. 
To calculate the correlation function for pairs of iden- 
tical particles we need to know the emission function 
S IW (x,p) describing the phase-space distribution of the 
particles emitted from the source. Since both the RWM 
and LEXUS provide explicitly only momentum space in- 
formation, we must try to reconstruct the corresponding 
coordinate space information from the description of the 
models given in Refs. [^J^]. From the emission function 
the correlator is calculated via [§-|lo| 



c iw {q,k) = i 
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where K = (pi + p2)/2 is the average pair momentum 
and q = pi — P2 the difference between the two observed 
on-shell momenta. Since pi and P2 are on-shell, q and K 
satisfy the constraint 
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The single-particle spectrum is calculated as 

E p ^-=P x (p) = J d 4 xS vw {x,p,E p ) . 
For the RWM it was given in the form Q 

J dYdpd<£f(p)d(Y L -\Y-Y \) 

V m± ( p-u 
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Here y,p±, 4> describe the momentum p of the measured 
particle, T is the temperature of the fireballs formed 
in the individual N-N collisions, and Y is their longi- 
tudinal rapidity which is distributed with a box distri- 
bution 0(Y L -\Y-Y o \) between Y -Y L and Y +Y L . u 
is the 4-velocity of the fireball, parametrized according 
to = (cosh p cosh Y", sinhpcos$, sinhpsin$, 

coshpsinhY") = 7(1,14) where p is the transverse rapid- 
ity of the fireball and $ its direction. This gives 

p ■ u=m± cosh(K — y) cosh p — p± cos($— <f>) sinh p . (5) 

The transverse fireball rapidity p is distributed with 



f{p) = f AB (p) 
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The width of this distribution depends on the collision 
system A+B via the random walk rule 



(N A + N B -2) S 2 



(7) 



where Na (N b ) is the average number of nucleons in 
nucleus A (B) hit by a nucleon from nucleus B (A). If the 
projectile is just one nucleon, Nb = 1. The parameter 8 
is obtained by fitting pA data (Nb = 1). 

If one tries to reconstruct S IW by comparing (||) with 
([I|) one is faced by certain ambiguities, although the 
appearance of a simple volume factor Vq in front of 
the integral over the fireball 4-velocity u in (Q) ex- 
cludes most functions S TVJ (x,p) with complicated x- 
dependences. The simplest assumption that the longi- 
tudinal and transverse fireball rapidities Y and p of the 
fireballs are not correlated with the fireball positions is, 
however, untenable. In this case the emission function 
factorizes, S rw (x,p) = F(x) ■ I(p) (with I(p) given by the 
r.h.s. of Eq. divided by Vo), and (|l|) yields a corre- 
lator which does not depend on K . Since all heavy-ion 
data show a clear dependence of the correlation function 
(in particular of its longitudinal width parameter Ri) on 
the pair momentum K |Ti]-[l| , such an emission function 
is excluded. 

More realistically one should at least implement the 
simple expectation that (in RWM language) fast fireballs 
live longer and fly farther before decaying than slow ones. 
This leads to a correlation between the decaying fireball's 
position x and velocity u, inducing also a correlation be- 
tween x and p via the Boltzmann factor and thereby a 
if-dependence of the correlator. We thus write 



5 rw (x,p)= J dYdpd<S>F(x;Y,p) 



x 703 cosh(F - y) exp ( - 
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with J d i xF(x;Y,p) = V 9(Y L -\Y-Y \). Azimuthal 
symmetry excludes a ^-dependence of F. We will 
use longitudinally boost-invariant space-time coordinates 
xf = (r,r],r,(p), withr 2 = t 2 ~z 2 , rj = \ \n[(t+z)/{t-z% 
and integration measure d A x — r dr dr\r dr dip. 

In the RWM |lj the transverse fireball rapidity is not 
correlated with its space-time position since the width 
Sab of f(p) is determined from the impact parameter av- 
eraged number of hit nucleons. Thus there are no trans- 
verse x-p correlations in the source. This is easily seen to 
lead to a correlator whose transverse width, given by the 
transverse HBT radius R±, does not depend on K . This 
contradicts the experimental observation of a clear and 
rather strong decrease of R± as a function of K± in 160 
A GeV Pb+Pb collision s [p" l| which cannot be explained 
by resonance decays [jfj 15 [ 



One can try to include the missing transverse x-p cor- 
relations by making the broadening parameter Sab r- 
dependent. Instead of using the impact parameter aver- 
aged number of hit nucleons from the original formulation 
of the RWM H|, 



N A = (2/3) ixrl 2R A n . 



(9) 
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(where ro — 0.8 fm is the nucleon radius, Ra = 1.12 A 1 / 3 
fm is the nuclear radius, and no = 0.17 fm -3 is the stan- 
dard nuclear density), we can take the actual value at 
distance r from the collision axis: 



N A (r) = (nr 2 )2y/R 



r 2 uq . 



(10) 



Via Eq. (Q) this yields an r-dependent width Sab and, 
via (^J), an r-dependent p-distribution Jab (fi) f). We will 
now study whether the resulting transverse x-p correla- 
tions can cause the measured -ftfj_-dependence of R±. 

Since the RWM does not take into account longitu- 
dinal energy degradation by multiple initial state scat- 
tering, we can factorize the p and Y dependence of the 
function F in (§): F(x; Y, p) = F'(x; Y) f(p, r). There 
is no correlation between the transverse fireball rapidity 
and its longitudinal position. Furthermore, without lon- 
gitudinal momentum loss the maximum and minimum 
longitudinal fireball rapidities Y~ ± Yl are not related 
to the actual number of hit nucleons and thus indepen- 
dent of r. Hence F'(x; Y) can be further factored into 
H(t], t; Y) ■ G(r, ip). The transverse HBT radius R± can 
be calculated from 11171 
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r sin 
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Since the radial integral couples G and /, we must know 
the functional form of G(r, <p) . We take G proportional to 
the number of collisions at distance r from the collision 
axis; for nuclear collisions at zero impact parameter this 
is given by 



G(r,ip) = G(r)<xy/R 2 A -r 2 y/. 



1 R 2 B — r 2 . 



(12) 



In the following we consider only symmetric collisions, 
Ra — Rb- The generalization to the case Ra ^ Rb is 
obvious. 

The evaluation of ([ll]) is complicated by the fact that 
Eq. (||) couples p to Y such that the function H(t),t; Y) 
must also be known. To obtain an as much as pos- 
sible model independent estimate for the behaviour of 
R± we use the following simplification: We go to the 
LCMS (Longitudinally Co-Moving System, y = 0) and 
fix Y = 0, thus taking into account pion production 
only from fireballs which rest in the LCMS. Since con- 
tributions from other fireballs are suppressed by a factor 
exp[— Mj_(coshY — 1) cosh p/T], we expect our simplifi- 
cation to provide a good estimate for the correct R±. 
The advantage of this approximation is that now the in- 
tegration in Eq. ( |Tl| ) factorizes in (r, 77) and (r, ip) . One 
finds 



J dp J r drr 2 C (r, p) 
JdpJrdr2C(r, p) 
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where 



C(r, p) = (R 2 A ~ r 2 ) (tt (4tt r 2 n Q ^R 2 
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x exp 
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( Mi \ (K x \ , x 
xexpl — — coshp j I I — smhp j . (14) 

Here the integrations over $ and ip have already been 
performed. Note that the r-integration in Eq. ([l3]) does 
not extend to Ra but only to Ra — £ where e is chosen 
such that at transverse distance r = Ra — £ there is 
still one nucleon-nucleon collision. This is ensured by 
requiring 5aa{Ra — s) — 0; if r > Ra — £ then Saa(t) is 
negative due to Eqs. (Q) and ([To|). 
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FIG. 1. The Mi -dependence of R± according to ([TJ]) for 
Pb+Pb (A=207) and S+S (A=32) collisions. 

For our computation we re-calculated S according to 
Eq. (|) (see ||). Using S pW = 0.3 § we got 8 = 0.217. 
Results for Rx according to the approximation (H|) for 
S+S (.4=32) and Pb+Pb (A=207) reactions are shown in 
Fig. [y. It is immediately clear that the calculated values 
for Rx underpredict the experimental data (which, at 
small Mj. , give Rx ~ 6 fm for Pb+Pb [[llj and « 4 fm for 
S+S |T^,[l3|]) considerably. It has been argued in |18| that 
the present data on Rx from Pb+Pb collisions require a 
strong transverse expansion since the measured values 
of Rx signal a source much bigger than the original Pb 
nucleus. The Rx -values in Fig. ^, however, are consistent 
with a transversally non-expanding source. Furthermore, 
the M^-dependence of Rx induced by the r— dependence 
of 5ab is much too weak and cannot reproduce the strong 
A/i-dependence seen in the Pb+Pb data . 

The wrong magnitude and if-independence of Rx are, 
however, not the only problems of the RWM/LEXUS 
models. Following |19| we now argue that they are also 
unable to correctly describe the dependence of the lon- 
gitudinal radius parameter Ri on the size of the collision 
system. 
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In RWM and LEXUS the full source of secondary par- 
ticles consists of smaller sources generated in individual 
N-N collisions. Let s c (x,p;xo) be the emission function 
describing such a small source where xo stands for the 
point where the N-N collision occurs. Translation in- 
variance gives s c (x,p;xo) = s c (x — Xo,p), and the total 
emission function can be written as 



S(x,p)= / d 4 x a(x ) s c (x - x ,p) , 



(15) 



where <t{xq) is the distribution of collision points re- 
sulting from the collision geometry. Its normalization 
/ d 4 xo a(xo) is fixed by the number of N-N collisions. 

We are now interested in the qi -dependence of the cor- 
relator. Inserting jl5| ) into (|l|), setting q± = 0, and boost- 
ing into the LCMS (# = = K t ) we find 



C(qi,K±) 



qj_=Kl=0 



\J d A x d 4 xo e mz a(xo) s c {x — xq, K) 
| J d 4 x d 4 XQ o-(xq) s c (x — xq, K) | 2 



/ d 4 xe~ lq ' z s c (x,K) 



J d 4 xs c (x,K) 



J d 4 xe~ lq ' z a(x) 



J d 4 xa(x) 



(16) 
(17) 



where z = r sinhry. This expression is typical for LEXUS 
type models where <j(x) is determined by collision geom- 
etry only. However, an additional if-dependence of a 
would not affect the following argument. The correlator 
would still factorize, and cr(x) will still be non-zero only 
in the region where N-N collisions occur. 

The correlator ( |l7j ) is in fact a product of two corre- 
lators associated with the sources s c and a, respectively. 
The longitudinal radius parameter of the whole source 
extracted from a Gaussian fit to (17) is thus given by 



R? 



R-l.a J 



(18) 



where Ri, Sc (Ri,a) is the longitudinal radius parameter 
associated with the source function s c (a). Ri. Sc can be 
extracted from correlation measurements in N-N colli- 
sions. Based on data published in |20|-|22|] this value was 
estimated in jl9| to be Ri, Sc ~ 1-7 fm. 

When considering different collision systems clearly 
only Ri a changes. It can be calculated using the model- 
independent expression p3] 



two nuclei are replaced by cylinders with the same vol- 
ume and (z 2 ). One finds fl9|| 



Rl,a = 
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where 7 w 10 is the Lorentz contraction factor in the 
N-N CMS at CERN/SPS energies. The factor l/(207 2 ) 
makes Ri, a rise by about a factor 6 more slowly than 
required by the data; for Pb+Pb collisions, after adding 
the contribution Rf Sc ~ 3 fm 2 , Eqs. ( p^ , ptj| ) underpre- 
dict the data |[l| for Rj by a factor 10. Hydrodynamic 
parametrizations 23 24] of the source function, on the 
other hand, can easily accommodate this rise by adjust- 
ing the lifetime To of the collision fireball. For Pb+Pb 
collisions at CERN one finds r « 8 fm/c J2|. This 
parameter measures the length of time during which the 
secondary particles created in the nuclear collision rescat- 
ter among each other. In scenarios with approximately 
boost-invariant longitudinal expansion Ri is directly pro- 
portional to this time. Models without rescattering (like 
RWM and LEXUS) do not lead to such a collective lon- 
gitudinal expansion and thus cannot reproduce Ri. 

In summary we conclude that the RWM and LEXUS 
models, which try to describe the hadronic momentum 
spectra in A-B collisions in terms of a linear extrapola- 
tion of pA respectively pp spectra, cannot reproduce the 
magnitude of the transverse radius parameter R± and its 
dependence on the transverse momentum nor the depen- 
dence of the longitudinal radius parameter Ri on the size 
of the collision system. The first failure is due to the lack 
of transverse expansion and transverse space-momentum 
correlations in these models. The second failure is caused 
by the lack of rescattering which would result in longer 
lifetimes and larger values of Ri through collective longi- 
tudinal expansion. Simple improvements of these models 
which preserve their original spirit as a superposition of 
N-N collisions cannot remedy these failures. Extensive 
final state rescattering among the produced hadrons, as 
implemented in other Monte Carlo transport models and 
in hydrodynamical source parametrizations, is essential 
to reproduce the correlation data. 

This work was supported by DAAD (B.T.), DFG, 
BMBF, and GSI (U.H.). We thank Andrei Leonidov for 
clarifying discussions. 
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where (f(x)) a = f d 4 x f(x) a(x)/ J d 4 xa(x). 

The compilation in Fig. 1 of Ref. |l^] shows that 
in A-B collisions Ri tCr rises approximately linearly with 
(A — l) 1 / 3 + (B— I) 1 / 3 , with a slope which is close to 0.5. 
The geometric distributions o~(x) for the collision points 
in RWM or LEXUS cannot (via Eq. (|l9|)) reproduce such 
a strong rise p9| . This can be seen immediately by cal- 
culating the r.h.s of (|l9| ) in an approximation where the 
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